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ABSTRACT 


* 


The  relative  motion  between  an  earth-based 
radar  and  a  satellite  object  is  of  interest  for  radar¬ 
tracking  studies,  signature  simulation  studies,  and 
for  developing  inverse  scattering  techniques.  Here, 
the  motion  analysis  is  accomplished  with  elementary 
vector  algebra.  This  approach  results  in  consider¬ 
able  simplification  of  the  final  equations  and  relation¬ 


ships. 
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SECTION  I 


« 


INTRODUCTION 

The  relative  motion  between  an  earth-based  radar  and  a  satellite 
object  is  of  interest  for  radar -tracking  studies,  signature  simulation 
studies,  and  for  developing  inverse  scattering  techniques.  This  relative 
motion  can  be  decomposed  into  two  component  motions:  the  radar -tracking 
motion;  and  the  rotation  of  the  satellite  about  its  center  of  mass.  The 
radar -tracking  motion  is,  in  turn,  caused  by  the  relative  motion  between 
the  satellite  in  its  orbit  and  the  radar  site  fixed  on  the  surface  of  a  rotating 
earth. 

The  equations  and  relations  required  in  consideration  of  overall  motion 
can  be  considered  in  terms  of  products  of  rotation  matrices  which  transform 
from  one  coordinate  system  to  another.  *"*"  The  approach  here  is  to  examine 
each  of  the  component  motions  separately  and  to  form  the  final  combination 
with  elementary  vector  manipulations. 

The  purpose  of  this  report  is  to  study  the  gross  effects  of  the  relative 
motion  on  radar  observations,  and  consequently  only  first-order  orbit  per¬ 
turbations  are  considered.  The  advantages  of  considering  higher  order 
perturbations  are  not  great  enough  to  justify  the  complications  which  they 
introduce. 


*  J.  F.  A.  Ormsby,  Motion  Simulation  of  Ground  Observed  Earth  Satellites, 
The  MITRE  Corporation,  W-7289,  Bedford,  Mass.  ,  4  September  1964. 

t  J.  F.  A.  Ormsby  and  S.  H.  Bickel,  A  Generalized  Radar  Output  Simulation, 
The  MITRE  Corporation,  W-7346,  Bedford,  Mass.  ,  19  October  1964. 


1 


Both  stabilized  and  torque-free  motion  of  the  satellite  about  its  center 
of  mass  are  considered.  In  all  cases,  expressions  for  the  relative  orienta¬ 
tion  angles  between  the  body  and  an  earth-based  radar  are  developed. 
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SECTION  II 


ORBITAL  MOTION 

There  are  four  different  motions  of  a  satellite  in  orbit: 

(1)  the  motion  of  the  satellite  which  is  restricted  to  move  in  an 
elliptical  path  due  to  the  inverse  square  law  for  radial  force; 

(2)  the  slipping  motion  of  the  perigee  point  within  the  plane  of  this 
ellipse; 

(3)  the  twisting  motion  of  the  plane  regressing  westward;  and 

(4)  the  change  in  the  size  and  shape  of  the  ellipse  due  to  atmospheric 

drag. 

Effects  (2)  and  (3)  above  are  caused  by  the  bulge  of  the  earth  at  the 
equator.  Five  orbital  elements  (see  Figure  1)  are  necessary  to  define  these 
motions: 

(1)  the  perigee  height  h,  which  is  the  distance  between  the  center  of 
the  earth  and  the  closest  point  on  the  ellipse. 

(2)  the  inclination  angle  6,  which  is  the  angle  that  the  orbital  plane 
makes  with  the  equator; 

(3)  the  right  ascension  of  the  orbit,  which  takes  the  form 

a  =  of  +  a;  t ,  (1) 

where  o-  is  the  right  ascension  at  time,  t  =  0,  and  is  the  rate  of  change 
in  the  westward  direction  of  the  right  ascension  due  to  the  earth  bulge; 

(4)  the  time  of  the  ascending  node,  which  takes  the  form 

2  3 

t  -  t  +  Tn  +  Dtn  +  D^n  ,  (2) 
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where  n  is  the  revolution  number,  T  is  the  satellite  period,  and  and 
D.}  are  drag  coefficients;  and 

(5)  the  argument  of  perigee,  which  takes  the  form 

<P  =  <PQ  +  k2  (f  '  sin2  <5)t>  (3) 

where  k  is  a  constant  which  specifies  the  rate  of  change  in  the  argument 

2 

of  perigee  due  to  the  earth's  bulge.  At  sin  6  =  4/5,  or  at  an  inclination  of 
approximately  63. 4  degrees,  this  rate  of  change  is  zero.  At  this  inclination, 
a  crossover  exists  where  the  rate  of  change  shifts  from  positive  to  negative. 

The  inverse  square  law  for  the  radial  force  field  from  the  center  of 
the  earth  results  in  an  elliptical  orbit  for  the  satellite  object  with  the  earth's 

A  A  A 

center  at  one  focal  point.  Let  e  f  e  ,  e  be  an  orthnormal  set  of  vectors  in 

1  u  O 


4 


4 


an  inertial  coordinate  system  with  the  origin  at  the  center  of  the  earth.  Let 
e  point  in  the  direction  of  vernal  equinox,  e  in  the  direction  of  the  north 

1  u 

pole,  and  e  in  the  direction  of  e  X  e  .  In  this  coordinate  system,  the 

u  u  1 

vector  equation  for  an  ellipse  becomes 

0  =  2  e^  £(a  •  a.  )  (cosE-e)  +  (  b  •  e.)  sinfij  ,  (4) 

i  =  1 

b  have  the 
;  ellipse, 
eccentric 

(5) 

where  t  is  the  time  of  perigee  crossing,  d  is  the  gravitational  constant, 
and  a  is  the  length  of  a.  Substituting  for  d  in  Equation  (5),  the  major 
semi-axis  of  any  satellite  is  given  in  terms  of  its  period  by 

a  =  205.  82T2//3  ,  (6) 

where  T  is  the  period  in  minutes.  The  ellipticity  of  the  orbit  is  given  in 
terms  of  the  semi-major  axis,  and  the  perigee  height  by 

e  =  1  -  h/a. 

The  semi-minor  axis  is  given  by 

,  / .  2 \  1/2 
b  =  a[  1-e  ) 


*  See  page  1,  first  reference. 


(?) 

(8) 


where  e  is  the  eccentricity  of  the  orbit,  and  the  vectors  a  and 
length  of  point  along  the  semi-major  and  semi-minor  axes  of  the 
respectively.  (Note:  a  points  in  the  direction  of  perigee. )  The 
anomaly.  *  which  is  related  to  time  in  Kepler's  equation,  is  E. 

/  a3  \  1/2  T 

t  "  t  =  (  -7  )  (E-e  sinE)  =  —  (E-e  sinE), 
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From  (4)  it  follows  that  the  length  of  the  vector  0  is  given  by 

1 0 1  =  a  1 1-e  cosE  | .  (9) 

Figure  2  illustrates  the  geometrical  relationships  between  these  orbital 
elements. 

In  order  to  complete  the  description  of  the  orbit,  the  vector  dot 
products  indicated  by  (4)  must  be  represented  in  terms  of  the  inclination 
angle  and  the  argument  of  perigee.  Let  a  be  the  angle  between  ascending 

A  A 

node  A  and  the  unit  vector  e  (see  Figures  1  and  3). 


i 

< 


Figure  2.  Basic  Ellipse  Orbit 
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£s  POLARIS 


Figure  3.  Earth-Orbital  Orientation 

From  the  geometry  of  Figure  3,  it  follows  that: 

(a  -  e^)  =  cos  a  cos  (p  -  sin  <p  sin  a  cos  6; 

(a*  e  )  =  sin  a  cos  <p  +  sin  <p  cos  a  cos  6; 

(a*  e  )  =  sin  6  sin  (p. 

u 

A  A 

Since  b  is  perpendicular  to  a,  it  is  only  necessary  to  replace  (p  by 
in  order  to  develop  the  corresponding  equations  for  b . 


(10) 

(U) 

(12) 

+  tt/2 
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SECTION  III 


SITE  ORIENTATION  AND  RADAR-TRACKING  MOTION 


Let  A  be  the  latitude  of  the  site  and  a  the  right  ascension  of  the 
s  s 

site  (i.  e.  ,  the  angle  between  the  site  longitude  and  the  vernal  equinox. 
This  angle  is  time  varying  due  to  the  earth's  rotation).  The  unit  vectors 
in  the  direction  of  the  zenith,  east,  and  north  as  measured  at  the  site 
are  given  by: 


Z  =  cos  X 

s 

(e„cosa:  +  e^sina  'j  +  sinA  en  ; 

VI  s  2  s/  s  3 

(13) 

E  =  -e„  since  +  e^cosce  ; 

1  s  2  s 

(14) 

N  =  -sinA  ( 
s  ' 

'  e„  since  +  e„since  )  +  cosA  e„ 
k  1  s  2  s/  s  3 

(15) 

(see  Figure  4  for  a  description  of  the  geometry  involved). 


Figure  4.  Site-Earth  Geometry 
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The  vector  from  the  site  to  the  object  is  given  by 


4 


4 


r  =  0  - 


roZ- 


(16) 


where  r  is  the  radius  of  the  earth.  From  Equation  (16),  it  follows  that 
e 

the  distance  or  radar  range  from  the  site  to  the  object  is 


r  = 


-  2r  (0  •  Z)  +  r 
e 


2 

e ' 


(17) 


The  evaluation  and  azimuth  angles  of  the  radar  (i.  e.  ,  the  angle  which  the 
target  makes  with  the  horizon  measured  at  the  site  and  the  angle  between 
due  north  and  the  target  measured  in  a  clockwise  direction  when  looking 
down  on  the  earth)  are  given  by: 


Tan  r  = 


(0-  Z)  -  r 
[|  0  I  “  -[(0-  Z)2]1/2 


(18) 


Tan  p  =  “4r-  >  0  <  p  <  2  7r  .  (19) 

0  •  N 

A  simple  and  useful  acquisition  scheme  would  be  to  test  for  the  zero 
elevation  angle.  Positive  angles  indicate  that  the  satellite  is  above  the 
horizon  and  visible,  while  the  satellite  is  below  the  horizon  for  negative 
elevations.  From  Equation  (4),  0-  Z  can  be  written  as 


0>  Z  =  OZ  cos  (E  -  EZ)  -  e(a-  Z), 


(20) 


where 


OZ 


=  \4a 


(a-  Z)2  +  (b-  Z)2 


Tan  EZ  =  — -| 
a .  Z 


(21) 
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Thus,  it  follows  from  setting  the  numerator  of  (18)  equal  to  zero  that  the 
horizon  plane  (i.  e.  ,  the  plane  tangent  to  the  earth  at  the  site)  intersects 
the  orbit  when  E  is  given  by  either 


E1  =  EZ  +  EC 


or  by 


E  =  EZ  +  2? r  -  EC, 


where 


EC  =  cos 


.  /  r  +e(a  •  Z) 
-1  /  e  '  ' 


OZ 


(22) 


(23) 


(24) 


The  first  solution,  E^  ,  indicates  the  last  point  where  the  satellite 
is  visible  before  it  sinks  below  the  horizon,  while  E  indicates  the  first 

Lt 

point  where  the  satellite  becomes  visible  above  the  horizon  (see  Figure  5). 

The  dot  product  between  0  and  any  general  vector  V  is  given  by 
(20)  and  (21),  where  Z  is  now  replaced  by  V.  Hence,  in  order  to  complete 
the  description  of  range,  elevation,  and  azimuth,  it  is  necessary  to  find 
the  vector  dot  products  between  the  set  Z,  N,  E,  and  the  orbit  axes  a  and 
b.  After  combining  (10),  (11),  and  (12)  with  (13),  (14),  and  (15),  the 
desired  products  for  a  become 


a*  Z  =  C,cosX  cos  a  +  a  sinX  sin  6  sin  w. 
Is  s 

a  •  E  =  -  C  sin  a , 

_  A 

a*  N  =  -  C.  sinX  cos  a  +  cosX  sin  6  sin  cp. 
Is  s 


(25) 

(26) 
(27) 


ft 
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where 


a 


=  (u)  +  w  t  +  «„  -  a  -  Tan  *  ( cos  6  tan  cp)  ; 
V  e  n/  1  n  V  > 


a >  =  angular  velocity  of  the  earth; 

oj  =  rate  of  change  of  the  right  ascension  of  the 
n 

ascending  node  due  to  earth's  bulge; 


=  right  ascension  of  the  site  at  t  =  0  ; 

a  =  right  ascension  of  the  ascending  node  at  t  =  0  ; 
n 


6  =  inclination  of  the  orbital  plane  ; 
cp  =  argument  of  perigee  ; 


X  =  latitude  of  site;  and 
s 


V  2  2  2 

cos  q>  +  sin  ip  cos  6. 


Figure  5.  Region  of  Satellite  Visibility 
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The  corresponding  equations  for  b  can  be  developed  by  replacing  cp 

by  cp  +  7r  / 2  in  the  above  equations. 

In  order  to  relate  the  radar  orientation  to  the  target  orientation,  it 

is  convenient  to  relate  the  radar  to  the  e  ,  e  ,  e  coordinate  system  or, 

more  specifically,  to  specify  the  latitude  and  right  ascension  of  unit 

A 

vectors  in  the  radar  range  direction  r,  the  radar  horizontal  H,  and  the 
radar  vertical  V.  From  inspection  of  Figure  4,  the  unit  vector  in  the 
radar  range  direction  is  given  by 


r  =  sin  t  z  +  cos  r  (cos  p  N  +  sin  p  E) .  (28) 

From  Equations  (13),  (14),  and  (15),  it  follows  that  the  latitude  of  r  is 
given  by 


sin 


X  = 
r 


r  •  e. 


sin  t  sin  X  +  cos  t  cos 
s 


X  cos  p 
s 


(29) 


and  the  right  ascension  by 


Tan  a 

r 


sin  t  sin  a  cos  X  +  cos  t  (sin  p  cos  a  -  sin  X  cospsina  \ 
S  S  V  s  s  s/ 

sin  t  cos  a  cos  X  -  cos  r  (sin  p  sin  a  +  sin  X  cospcosa  \ 
s  s  V  s  s  s/ 


(30) 


The  latitude  and  right  ascension  of  the  radar  vertical  vector  can  be 
found  by  replacing  r  by  t  +  ?r/2  in  (29)  and  (30),  while  the  horizontal 
vector  is  defined  by  letting  t  =  0  and  p=p+7r/2.  In  this  case,  the 

A  A  ^ 

vectors  r,  V,  and  H  form  a  right-handed  system. 
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SECTION  IV 


TARGET  MOTION 


» 


* 


In  general,  three  Eulerian  angles  are  necessary  to  specify  the 
orientation  of  an  earth-based  radar  with  respect  to  the  satellite.  They 
are: 

(1)  aspect  angle,  the  angle  between  an  axis  fixed  in  the  satellite 
(satellite  axis)  and  the  radar  line-of-sight; 

(2)  polarization  angle,  the  angle  which  the  projection  of  the 
satellite  axis  on  the  plane  normal  to  the  radar  line-of-sight 
makes  with  the  radar  horizontal  polarization  axis;  and 

(3)  roll  angle,  the  angle  which  denotes  rotation  about  the  satellite 
axis. 

These  angles  are  labeled  9,  (3,  and  £,  respectively,  in  Figure  6. 


Figure  6.  Euler  Angle  Geometry 
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B  points  in  the  direction  of  the  satellite  axis,  and  £  is  referenced  from  a 
unit  vector  in  the  r  X  B  direction  to  a  fixed  axis  B  in  the  body  perpen- 

X 

dicular  to  B.  The  horizontal  axis  H  is  taken  in  the  r  X  V  direction 

A 

(i.  e.  ,  H  is  pointing  east  when  the  radar  is  pointing  north).  It  follows 
from  the  geometry  of  Figure  6  that  these  angles  are  given  by: 


cos  9 

=  r  • 

B, 

0  <  0  <  7r  ; 

(31) 

Tan  /3 

A 

B- 

B- 

A 

V 

*  j 

H 

0  <  /3  <  2 it  ; 

(32) 

B  .  r 

Tan  £  =  z - 1  ,  0  <  £  <  2t r  .  (33) 

B  X  B  .  r 

In  order  to  work  out  these  relationships  explicitly,  the  satellite 
motion  about  its  center  of  gravity  must  be  studied.  The  principal  types  of 
stabilized  and  unstabilized  motion  are  summarized  herein. 

The  six  principal  types  of  satellite  stabilization  are: 

(1)  spin,  the  satellite  is  fixed  in  inertial  space; 

(2)  earth  center,  one  axis  of  the  satellite  points  toward  the  center 
of  the  earth; 

(3)  earth  horizon,  the  satellite  is  stabilized  with  respect  to  the 
local  horizon; 

(4)  magnetic,  an  axis  of  the  satellite  points  along  the  lines  of 
magnetic  force  (often  the  magnetic  field  is  used  for  damping 
for  earth  center  stabilization);  and 


4 


4 
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(5)  inertial  guidance,  the  satellite  orientation  is  governed  by  some 
predetermined  program  or  by  signals  from  ground-based 
transmitters. 

The  four  principal  types  of  unstabilized  motion  are: 

(1)  tumbling, 

(2)  spin, 

(3)  precession,  and 

(4)  precession  with  nutation. 

By  combining  Equations  (31)  through  (33)  with  the  body  motion 
equations  for  a  particular  satellite,  the  three  orientation  angles  (aspect, 
polarization,  roll)  can  be  found  explicitly.  For  example,  consider  an 
earth-centered  stabilized  body.  In  this  case,  the  vector  B  points  along 
0  and  consequently 


cos 


r  sinr 
e 


(34) 


(35) 


The  roll  angle  will  depend  upon  the  type  of  stabilization  assumed  for  B  . 

X 

One  possibility  for  earth  center  stabilization  would  be  to  take  in  the 
direction  to  the  unit  vector  normal  to  the  orbital  plane.  In  this  case,  from 
Equation  (33) 


n  *  Z 

Tan  k  =  i  ~ — r-r- - -  -  -  —  ’  (36) 

(0  •  a)  (b  •  Z)  -  (0  •  b)  (a  •  Z) 

A  A 

where  n  =  a  X  b.  Explicit  formulas  for  the  dot  products  indicated  in  (36) 
were  presented  earlier. 

Another  common  method  of  stabilization  is  spin-stabilization.  Here, 
the  B  axis  is  considered  as  fixed  in  space.  Taking  A  ^  and  a  as  the 
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latitude  and  the  right  ascension  of  the  body  axis,  respectively  (which 
corresponds  to  the  angular  momentum  vector  for  spin- stabilized  bodies), 
the  aspect  angle  is  given  by 


cos  0  =  sin  X  sin  X  +  cos  X  cos  (a  -at  ) .  (37) 

pr  p  r  p  '  r  p'  ' 


Hence,  X^  and  q;^  correspond  to  the  latitude  and  right  ascension  of  the 
radar  in  inertial  coordinates.  *  Similarly,  the  dot  products  in  (32)  for  the 
angle  /3  may  be  determined  from  (37)  by  replacing  X^  and  a  by  the 
corresponding  latitudes  and  right  ascension  for  the  radar  horizontal  and 
vertical  vectors. 

Since  the  angular  momentum  vector  is  fixed  in  space  for  torque  - 
free  motion,  the  0  given  by  (37)  is  the  angle  between  the  radar  and  the 
angular  momentum  vector.  Torque-free  motion  of  an  axially  symmetric 
body  results  in  precession  about  this  vector.  If  we  take  0^  and  0  as  the 
precession  cone  angle  and  precession  rate,  respectively,  then  the  radar 
aspect  angle  becomes 

cos  0  =  B •  r  =  cos  0  cos  0  +  sin  0  sin  0  cos  (d>t+  <p  -0,')  ,  (38) 

pr  p  pr  p  V  o  1/  ' 


where 

Tan 


sin  X  cos  X  -cos  X  sin  X  cos  (a  -a  \ 
_ r  p _ r _ p  V  r  p/ 

cos  X  sin  ( a  -a  X 
y  v  r  p/ 


(39) 


and  0  is  the  initial  location  at  time  t  =  0  of  the  body  axis  on  the  preces- 
o 

sion  cone.  The  particular  cases  of  spin-stabilized  or  tumbling  motion  may 
be  obtained  by  setting  the  precession  angle  equal  to  0  or  7r/2  ,  respectively. 


*See  Equations  (29)  and  (30). 
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(5)  inertial  guidance,  the  satellite  orientation  is  governed  by  some 
predetermined  program  or  by  signals  from  ground-based 
transmitters. 

The  four  principal  types  of  unstabilized  motion  are: 

(1)  tumbling, 

(2)  spin, 

(3)  precession,  and 

(4)  precession  with  nutation. 

By  combining  Equations  (31)  through  (33)  with  the  body  motion 
equations  for  a  particular  satellite,  the  three  orientation  angles  (aspect, 
polarization,  roll)  can  be  found  explicitly.  For  example,  consider  an 
earth-centered  stabilized  body.  In  this  case,  the  vector  B  points  along 
0  and  consequently 

cos  0  =  ~  ( | 0 | -  r^  sin  r  )  (34) 

P  =  \  •  (35) 

The  roll  angle  will  depend  upon  the  type  of  stabilization  assumed  for  B  . 

X 

A 

One  possibility  for  earth  center  stabilization  would  be  to  take  B^  in  the 
direction  to  the  unit  vector  normal  to  the  orbital  plane.  In  this  case,  from 
Equation  (33) 

_  .  n  •  Z 

Tan  *  =  — ; — — - _  7  -  „  ’  (36) 

(0-  a)  (b-  Z)  -  (0-  b)  (a-  Z) 

A  ^  ^ 

where  n  -  a  X  b.  Explicit  formulas  for  the  dot  products  indicated  in  (36) 
were  presented  earlier. 

Another  common  method  of  stabilization  is  spin-stabilization.  Here, 

the  B  axis  is  considered  as  fixed  in  space.  Taking  X  and  a  as  the 

P  P 
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latitude  and  the  right  ascension  of  the  body  axis,  respectively  (which 
corresponds  to  the  angular  momentum  vector  for  spin-stabilized  bodies), 
the  aspect  angle  is  given  by 


cos  9  =  sin  X  sin  X  +  cos  X  cos  (a  -a  ) . 

pr  p  r  p  '  r  p' 


(37) 


Hence,  X  and  a  correspond  to  the  latitude  and  right  ascension  of  the 
r  r 

radar  in  inertial  coordinates.  *  Similarly,  the  dot  products  in  (32)  for  the 

angle  f3  may  be  determined  from  (37)  by  replacing  X  and  a  by  the 

r  r 

corresponding  latitudes  and  right  ascension  for  the  radar  horizontal  and 
vertical  vectors. 

Since  the  angular  momentum  vector  is  fixed  in  space  for  torque - 

free  motion,  the  9  given  by  (37)  is  the  angle  between  the  radar  and  the 
pr 

angular  momentum  vector.  Torque -free  motion  of  an  axially  symmetric 
body  results  in  precession  about  this  vector.  If  we  take  0^  and  6  as  the 
precession  cone  angle  and  precession  rate,  respectively,  then  the  radar 
aspect  angle  becomes 

cos  9  =  B*  r  =  cos  9  cos  9  +  sin  9  sin  9  cos  (<bt  +  <p  -<£,')  ,  (38) 

pr  p  pr  p  V  o  ^1/  ’  '  ’ 


where 
Tan  <p  i 


sin  X  cos  X  -cos  X  sin  X  cos  fa  -a  ) 
_ r _ p _ _r _ p  V  r  p/ 

cos  X  sin  (a  -a  X 
y  v  r  p' 


(39) 


and  4>°  is  the  initial  location  at  time  t  =  0  of  the  body  axis  on  the  preces¬ 
sion  cone.  The  particular  cases  of  spin-stabilized  or  tumbling  motion  may 
be  obtained  by  setting  the  precession  angle  equal  to  0  or  7r/2  ,  respectively. 


4 


V 


*See  Equations  (29)  and  (30). 
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